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Abstract 

We discuss unstable particle mixing in CP- violating weak decays. It is shown that 
for a completely degenerate system unstable particle mixing does not introduce 
a CP- violating partial rate difference, and that when the mixings are small only 
the off-diagonal mixings are relevant. Also, in the absence of mixing, unstable 
particle wave function renormalization does not introduce any additional effect. 
An illustrative example is given to heavy scalar decays with arbitrary mixing. 



1 Introduction 

The smallness of — Ks mass difference allows us to have an access to rare processes 
such as CP violation. Up to now the only established experimental evidence of CP violation 
comes from the mixing of the unstable particles and . 

Earlier studies of unstable particle mixing followed two physically equivalent paths. One 
is due to Weisskopf and Wigner 0, in which one introduces an effective complex mass matrix. 
The evolution of the system is determined by the standard time-dependent Hamiltonian 
formalism 0. The other is due to Sachs 0], in which one studies the dynamics of the 
complex pole of the kaon field propagator. The Hamiltonian method is expressed directly in 
terms of the measured quantities and is therefore more transparent from a phenomenological 
viewpoint. On the other hand, the propagator method arises naturally in the context of 
quantum field theory, and hence is more easily adapted to fundamental gauge theories of weak 
interactions. Both approaches are phenomenological, having difficulties handling ultraviolet 
divergences arising from higher order corrections. In spite of these fundamental difficulties, 
the phenomenological formalisms have been very successful. They provide the standard 
descriptions for the study of unstable particle mixing. 

The advent of renormalizable gauge theory provides a connection between the parameters 
of a phenomenological formalism and the parameters of a given fundamental theory. In this 
paper we would like to study these connections for unstable particle mixing in some detail, 
focusing on CP-violating processes. We will adopt an approach that combines the two 
methods mentioned above. 

Instead of introducing a complete renormalization prescription, our immediate goal is 
more modest. In the next section we discuss some general properties of S*- matrix elements 
in the presence of unstable particle mixing. The results of this analysis turn out to be very 
useful for simplifying Feynman diagram calculations. 

In section 3 we study the relationship between the unstable particle mixing and antiparti- 
cle mixing. For simplicity, we only focus on scalars. A simple formula valid for small mixings 
is derived for CP-violating partial rate differences. 

The formalisms developed in section 2 and section 3 are applied to a simple example 
of baryogenesis by heavy scalar decay. The results are shown to agree with the published 
results obtained directly from Feynman diagram calculations. This part is presented in 
section 4, followed by a discussion in section 5 of large mixing and renormalization. Our 
conclusion is presented in section 6. We give two appendices to present some technical 
details: one discusses the renormalization of unstable particle mixing and the other shows 
how to diagonalize an arbitrary n x n complex matrix. 



2 General formalism 

Consider the weak decay of a set of unstable particles (pa produced at t = 0, where the index 
a = 1,2, ... labels different flavor of 0. Suppose the lowest order amplitude of 0^ decaying 
into a final state \Ff) is given by 
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By CPT invariance, to first order of -ffweak a replacement of Ff and 0^ by their antiparticles 
Ff and 4>a corresponds to change Tfa to T^^. By the superposition principle, the weak 
amplitudes at a later time t are 

Tfait) = J2TfbV,-'V,ae-'^^\ (2) 

b,c 

Tfait) = Y.T}bKc'Vcae-'^^\ (3) 

b,c 

where V and V are the mixing matrices 

l0a) = Ka|0:.), (4) 
\^a) = V,aWc). (5) 

and 10^), 10^) are the eigenstates of an effective Hamiltonian, which is correct up to second 
order in H^eak- |0c) 10'^) will be referred to as the eigenstate of propagation hereafter. 
The eigenvalue of |0^) and |0^) is Uc = nic — i'yc/'^ in the rest frame 0[,, where rric and 7c may 
be interpreted, respectively, as the mass and width of 0^ and 0'^. 

If \Ff) {\Ff)) belongs to the unstable particle set, then Tff, and T/b zero unless 
f = b, and Eqs. (Q) and are useful for the study of time distribution of CP asymmetry 



in oscillation |Tl], |12|. Except for situations in which the particles are stable or the mass 
and decay matrices commute, V and V are in general not unitary. If both CPT and CP 
are conserved, V = V and Eqs. (0) and (|^) are the natural generalizations of the known 



formalism |T3[ describing the so called 'mix-and-decay' phenomena. The formalism of |jl3 



has been employed in the study of unstable neutrino oscillation [T^, |T5[ . For practical purpose 
the mixing matrix has so far been approximated as unitary. Such an approximation is not 
always justifiable. 

Our main interest is in the unstable particle mixing effect in the time-integrated rate 
difference 

A;, = r(0„-.F^)-r(0,-.F;). (6) 

It should be pointed out that, in addition to the mixing, the rate difference may also depends 
on final-state interactions, and it is not always possible to separate them if these two effects 
are comparable. A discussion on particle mixing is always warranted, however, unless its 
effect is negligible. For recent discussions on the final-state interaction effect see Refs. 0- 

In order for A fa be nonzero it is necessary that CP be violated and to have significant 
non-trivial CP-conserving phases. According to Eqs. and (|^), a CP-conserving phase 
may arise either from the mixing matrix V or V or/and from the evolution phase e"*'^''*. For 
a completely degenerate system, i.e., Ua = uj, the evolution phase factors 



Tfait) = e~^-%a, (7) 
Tfait) = e— *T;,. (8) 



Although a nonreal Tfa implies CP violation, the contribution to the rate difference from 
mixing is seen to vanish because |Tya(t)p = \Tfait)\^. 



3 



This result has an important imphcation in searching for mechanisms for baryogenesis 
16| . It has previously been suggested that a degenerate unstable system might provide a 



resonance enhancement to the CP- violating partial rate difference from mixing. Our analysis 
shows that the outcome would be the opposite if the degeneracy is complete. Instead of a 
resonance enhancement, we expect that contributions to the generation of baryon number 
asymmetry from unstable particle mixing is highly suppressed whenever the particles in 
question have nearly equal mass and lifetime, and is zero in the complete degenerate limit. 

One special example of a completely degenerate system is that the set contains only 
one particle, i.e., a = b = c = 1. In that case (0) and (|^) show that particle instability 
itself does not contribute to the CP-violating partial rate difference. The same statement 
applies to systems with an arbitrary number of unmixed particles. In terms of perturbation 
theory these results imply that, in the absence of mixing, unstable particle wave function 
renormalization does not affect A fa- 

Had we considered the rate difference of the eigenstates of propagation, i.e., Aj^ = 



~^ Ff) — T{(j)'i — > Ff), we would have reached a different conclusion [|1^]. This is 
evident by dropping the matrix V (not V~^) from and V from (^). Phenomenologically, 
neglecting V and V in (H) and (|^) corresponds to ignoring mixing in particle production. 
As pointed out earlier ^ , the propagation eigenstates cannot be regarded as physical in 
the sense that they cannot be directly produced nor detected pO[ . As a result 
physical observable. 

Another important feature of Eqs. (0) and (|^) is that, when the mixings are small, the 
diagonal mixings and phases are irrelevant for A/q. Indeed, for small mixings, we expand V 
and V as 



Aj„ is not a 



4a + AKa 
5ca + AV^ca 



(9) 
(10) 



where the diagonal phases Uc and are real. The elements in AV and AV are assumed to 
be much smaller than unity. To first order of AV and AV we have 



b 

ffait) = e-^-*T;, + ^T;,AH, 



-iuj}jt 



— e 



-iuJat 



(11) 

(12) 



The second terms in Eqs. (^) and ( [T^ ) vanish for diagonal mixings. Therefore, neglecting 
off-diagonal mixings the CP-conserving phase e"*'^"* factors and |Tja(t)p — |Tya(t)P = 0. 
Hence, for small mixings the diagonal phases a^, cxa and Uat do not enter into the determi- 
nation of Afai only the off-diagonal mixings, i.e., h ^ a, are relevant. This is similar to a 
result of Wolfenstein's for final-state interactions . Thus, for the calculation of A ja one 
does not need to consider flavor- conserving one-particle-reducible diagrams. 



3 Relations between particle and antiparticle mixings 

CPT invariance provides a relationship between V and V . In quantum mechanics this would 
be obtained by studying the Hamiltonian. The existence of the 'standard model' underlines 
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the usefulness of a field theoretical analysis. In field theory the relation between V and V 
can be easily obtained from particle propagator. Consider situations in which 0^ is a scalar. 
At tree level, the propagator of 0^ is 



p2 



m 



(0)2 



5, 



(13) 



where m^") is the bare mass of (pa, and we have implicitly assumed that \(j)a) is an eigenstate 
of zeroth order of H^eak- Including one-loop corrections the following changes occur in (|1^) 



m 



ba 



ba' 



where A^^ is the square of an effective complex mass matrix 



Ml 



Ml 



(14) 



(15) 



in which Ai is the effective mass matrix and = |(A4r + TAi), where T is the decay 



matrix, and F /4 has been neglected. Both /A and T are hermitian, but A4 is not 
For a given P^, Ai^ can be diagonalized by a transformation 



QM^Q-' 



ba 



{ml - imaJa)Sba, 



where Q is a complex matrix. The one-loop regularized propagator can be written as 

lA^^J (P2) = «Q,-; _ ^2 ^ ^1 -Ig^^^ 



(16) 



(17) 



where for simplicity we have neglected terms of order 7^. Evidently, i — + imcjc 
is the propagator of 0^. It follows that 



Vta- 

and thus the relationship between V and V is 

Vba 



Qbai 
Qab ' 



(18) 
(19) 



(20) 



For stable particles Ai^ is hermitian, Q is unitary and hence so are V and V, and (pO]) reduces 
to the known result V = V*. Substituting (pOf ) into (|) one sees that the time-dependent 



mixing matrix in the antiparticle decay is Y,c ^ac^^~^'^''^^cb, which differs from that in the 
particle decay (|^) by exchanging the indices a and b (a consequence of time reversal). 

It is important that Ai^ is momentum dependent. This is necessary if the orthonormality 
conditions are to be maintained for both |0a) and |0'^) In practice this does not introduce 
any additional complication, as P^ is always fixed by the on-shell condition once the initial 
state is specified. 

We now focus on a case of special interest, small mixing. By small mixing we mean that 
(1) the width differences of the particles are much smaller than their mass differences and 
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(2) the off-diagonal elements in Ai^ can be treated as a perturbation. In that case AV and 
AV" defined by (|^) and (p!0D, may be separated into their dispersive and absorptive parts 



{A) 
ba 5 



(21) 
(22) 



AVha 

Phenomenologically, AV^^"^ and AV^f"^ correspond to the mixings arising from the mass and 
decay matrices, respectively. The hermiticity of the mass and decay matrices then implies 



AK 



(D,A) 



ba 



The solution to (pOD satisfying the constraints of (^31) and (|^) is 



This is a much simplified version of (^0]), valid for small mixings. 

With ( p5|) one can have a simple expression for A fa- From (|l^ 
that the time-differentiated CP-violating partial rate difference is 



(23) 
(24) 



(25) 



and f^Sf) we find 



\Tfa{t)\' - \ffa{t)\' = 4^p,{sin/36e2i--"* - sin[Re(o;, - ujb)t + (3b]e'^^'^'^+^^^'} , (26) 



where 



Pb = 
tan Ph 



Im 



+ 
/Im 



;im(AHl^)T;,r;, 



(27) 
(28) 



In the absence of CP violation AV^'^'^'"^) and Tja are real, Pb = and hence Eq. ( PB| ) vanishes. 
As pointed out earlier, Eq. (p6D shows explicitly that for diagonal mixing, i.e., b = a, the 
phase makes no contribution. If the mass and decay matrices commute, AV^*^"^-' = and 
reduces to 



|T/a(t)p-|T;,(t)p = -4^Im 



sm 



Re(a;a - ujb)t 



-,Iia{ula+OJb)t 



(29) 



The result given by (|29| ) can also be obtained from a general formalism developed recently 
by Gronau and Rosner for a 2 x 2 system. 

Within our approximation all oscillatory terms are integrated to zero. Hence 



A/a = [ dn r ^[\Tf^{t)\' - \ffM' 

J J TVi 



4 / rffi ^ Im 



(30) 



where l/r^ = —2lm.uJa is the width of the particle and / dQ represents a phase space sum. 
Eq. ( pO]) shows that for the calculation of A fa one only needs to consider off-diagonal mixings 
in the decay matrix. 
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For small mixings it is easy to show that 

E^/<^ = 0- (31) 

/ 

This relation follows from unitarity of S-matrix, which in the present situation implies 

Y: T},Tf,\p.=^. oc f L oc AVi^\ (32) 
/ 

Thus, S/.fe^a ^^Ih^'^afT^hf IS real and Eq. (0) follows as a consequence. Since the final-state 
interaction contribution to the total rate difference of a particle and its antiparticle is known 
to vanish 0, the result of (^) is in accordance with the usual expectation that a particle 
and its antiparticle have the same lifetime p3|. 



4 An illustrative example for small mixing 



In this section we show how to apply ( |30D to model calculations. Consider a system containing 
two heavy scalars Sa,a {a = 1,2) and a as a color index. The lowest order interaction of the 
system is given by 



(33) 



where u, d and e'^ are the charged fermions of the first generation, which are considered as 
massless. This interaction can arise from an SU{b) GUT theory with two 5-plets of Higgs 
fields. Since the two final states into which the heavy scalars can decay have different baryon 
number, (|33| ) provides an interaction for baryogenesis via heavy scalar decays. 

For convenience, we choose as a basis Sa^a the mass eigenstate fields and select = 
\dRu'j^ and IF2) = le/ju^). The lowest order transition matrix for IS'1^2) \Fi,2) is 

T = c{ ^^12 \ ^ ^34) 

Y G21 G22 / 

where C is an over all normalization constant. The factor \/2 in Tn and T12 is introduced 
to account for the effect of summing over the indices of e"^^ in calculating the squares of the 
elements. 

To determine AV^*^"^-' we consider corrections to the scalar propagators to second order of 
Lj. A simple calculation shows that the regularized elements in the complex-mass-matrix 
square ([T5|) are (remember and F^ are momentum dependent) 



Ml 



f2 

ba 



m 



167r2 



2Glf^Gia + G2hG2a 



2 3 
A-n ^ 2 



IE - In 



P2 



IGtt 



2Glf,Gia + G*2i^G2a 



(35) 
(36) 



where n is the dimension of regularization, /ig is the ultraviolet cut-off and 7^; is Euler's 
number. Following the standard technique 0] we find that the matrix which diagonalizes 
according to (^) is 



V 












(37) 
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where 6 = KeO + ilmO is complex 

tan 261 = 2- 



M2 + ^|f? 



121 



(38) 



with the real quantity = {Pyi6n^)\2Gl^Gi2 + ^21^22! [43;^ + | - 7i? - HPyATinl)]. 
The phase 5 is real and determined by 



6 = - are 
2 ^ 



(39) 



The eigenvalues are 



2 - «"^i,27i,2 



l{« + M22)-^(f?i + f: 



In the small mixing limit, i.e., 



M/i -M,^,l > |r 



221 



Ml 



one has from 



that 



11 -^2215 l^"12h 1-^ 12h 



-Im^e^^-^ 

l^Q^-^25 Q 



where 



Im^e 



125 



Im^e" 



i25 



(40) 



(41) 



(42) 



(43) 



Substituting (|34D and (^2|) into (|30|) and making use of the relation = 
which is valid to zeroth order in £/, we obtain for Si Fi^2, in which = ml, 



22 



m 



A 



11 



-A 



21 



0.1 ml 
271 {ml — ml 



-Im 



ll'-J12'-J2l'-^22 



1.25 



(44) 



where fli = mi/lGvr is a phase space factor for S'l. Also, for the S2 decays we have P^ = m\ 
and 



A 



12 



-A 



mo 



22 



2tt {ml — ml) 



Im 



ll{_r^2*-^2l'-^22 



(45) 



where Q2 = m2/lQTT. These results are in complete agreement with those obtained from a 
direct Feynman diagram calculation |^. Contributions to Aja from vertex corrections can 
be found in m, |2l 
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5 An illustrative example for large mixing 

In the small mixing limit (|4l|) one can simply use the regularized (but not renormalized) 
formalism to compute Aj^. The unphysical quantities in do not enter. However, to go 
beyond this limit we must introduce a renormalization prescription to remove the divergences. 
A pedagogic introduction illustrating how this may be done for the example discussed above 
is given in Appendix A. 

From a phenomeno logical viewpoint, the mixing phenomena under consideration is de- 
termined by the renormalized interaction lagrangian [EHI 



(46) 



where g fa [f = 1,2) are the renormalized couplings in the weak eigenstate basis. The square 
of a renormalized complex mass matrix is 



Ml; Ml, 



ill 12 

p2* 1^2 
-•^12 ■'-22 



(47) 



where M^^ and gfa are the parameters of the model determined experimentally (at some 
scale) . The parameters in F^^ are calculable 



f2 

ba 



p2 



167r2 



A convenient form for the matrix which diagonalizes ( ^T]) according to (|1^) is 

cos6' 



V 



sm I 



-i<f> 



cos 9 



(48) 



(49) 



Again, the mixing angle 6 is complex 



tan^ 29 = 4- 





'Ml, - iVl, 


'Ml;-iVll 




'{Ml, - Ml) 




^1 



The phase is also complex determined by 

-i2<j> 



Ml; - iTf, 



Ml, - iT^ 
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The eigenvalues are 



(50) 



(51) 



UJ 



1,2 



"^1,2 



:7i,2 



Ml 1 + M22 



22 



(ma - MI2) - <f ?i - fl,)] + a(mI, - if I,) [Ml; - ifli) . (52) 



If aigM? 



I, — argF^2) = 2aigMl, = 26 is real. Compared to the regularized formula 
one sees that Eq. (|^) differs from 



only by a diagonal phase matrix. This difference 
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has no physical significance. The diagonal phase matrix can be removed by a suitable choice 
of phase convention. 

Let us now turn back to Ajq. Eqs. ( ^81) to (|5^ ) enable us to compute Aja for arbitrary 
mixing angle and phase. The results are only limited by the validity of perturbation expan- 
sion of Cj. Here one should use the general formulas (0), and (^). In applying Eqs. 
(49) to (p2D one should also be very careful whenever the parameters are near the branch 



cuts in the complex parameter space. 

As an illustration we consider that the initial (renormalized) state is an eigenstate of 
mass matrix, i.e., ~^ 0; but not an eigenstate of the decay matrix. In that case is real 
and 



We find 



sin 29e 



-i(f> 



2* 



^11 



'^9119*12 + 9219*22 



(53) 



An 



-A 



21 



fiimpm 


9ll9l292l922 


271 


ul - ul 


2 



dt, 

Ti 



-Re 



)( 



)( 



,(54) 



where c = cos^^ and s = sin 9. Also, 

Ai2 = — A22 





9ll9l292l922 


27r 


2 2 


2 



T2 



-Re 



In ( |5^ ) and (0) the mixing angle 6 is arbitrary. 

In the limit in which the mixing angle is small due to |r^^ — 
Eqs. (B31) and (ESI) are simplified to 



An 
A12 



-A21 

-A22 



27T m{ 



ml 



27r ml 



2 



Im 



Im 



9ll9l292l922 
9ll9l292l922 



ml 



221 » wt^ 



mt. 



r2 _ r2 
11 J- 22 



11 22 



.(55) 



/if2 1 ir2 I 



(56) 
(57) 



Compared to (^^ and (|45|) ( after renormalization Gfa are replaced by (?/„), these results are 



further suppressed by a ratio {m\ • 



■m 



2\2 



/(r 



<^ 1. In particular, all the CP- violating 



partial rate differences vanish in the limit mi = m2. This is in contrast to what one might 
have expected, based upon a naive extrapolation from ( PD and (^S|). Since Ajq depends on 
the imaginary part of the mixing angle (|30|), it must vanish whenever Im6' = 0. For the case 
at hand, 6 is real when = and mi = m2- 

6 Conclusion 



We have discussed a formalism for unstable particle mixing base upon one-loop renormaliza- 
tion of field theory, with emphasize on its applications to CP-violating physical processes. 
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Among various interesting results, we have found that, for a completely degenerate system, 
unstable particle mixing does not contribute to the CP-violating partial rate difference. In 
particular, in the absence of mixing, unstable particle wave function renormalization does 
not introduce any additional effect for the CP asymmetry. When the mixing is small we 
show that only the off-diagonal mixings and phases are relevant for CP violation. 

We have used a simple example to show how to apply the formalisms developed in this 
paper in model calculations for arbitrary mixing. The basic steps for renormalizing unstable 
particle mixing are also outlined. 
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Appendix A 

It is most convenient to introduce a renormalization prescription in the weak eigenstate 
basis. We parameterize the interaction lagrangian as 



{51 



where gf} (/ = 1, 2) are the weak eigenstate bare couphngs. The one-loop regularized inverse 
propagator can then be written as 



+ 



1 . . ir„2 



167r2 



^g^gia + g*2b92a P\h - l - In ^ + /i^ + tnP'eiP') , (59) 



where m^^J'^ is the square of the bare mass matrix, /q = + f — 7s + ln47r, gja (/ = 

1,2) is the renormalized couplings (to the lowest order gfa = gf}) and 6 is the 6'-function. 
S^(P^,/i^) is the renormalized self-energy 



'^glbdla + g*2b92a 



167r2 

which has the standard features 



Sf„(/i2,/i2) = 0, 

^s£(p^/.2)lp.=,. = o, 



(60) 

(61) 
(62) 



and yU^ is the subtraction point. In practice it is convenient to choose /i^ be the invariant 
mass square of the initial particle. 

Following the standard method the renormalized inverse propagator is 



where r^^(P^) is given by ( ^81) and 



(63) 



(0)2 



ba 



167r2 



'^glbdla + g2b92a 



^— J2{ [2gl,g,, + g;,g2c] + Ml [2gl^gia + g*2,g2a] } [/. 



167r^ 



1-ln^ 



(64) 



is the square of the renormalized mass matrix. The renormalized weak eigenstate fields are 
related to the bare fields by 



167r2 



'^glbdla + glb92^ [/o - 1 - In —1 ]Sb,a- 



(65) 



The divergencies in the interaction lagrangian associated with wave function renormaliza- 
tion and vertex corrections are finally absorbed into coupling constant renormalization, i.e., 

(0) 

9 fa 9fa- 
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Once the renormalized parameters are determined at a given scale /i^ by experiments, 
the separation of the divergent term /o from and /i^ in (|60D, (|6^ ) and (^) allows us to 
predict these parameters at any other scale. 

Appendix B 

In this appendix we outline the basic ideas of diagonalizing an arbitrary n x n complex 
matrix Ti. Ti can always be diagonalized by a bi-unitary transformation, i.e., VLHVji ^ = 
Afo^ba, where Vl and Vr are unitary and the eigenvalues Aq are real. However, this procedure 
is not useful for us unless V^V^ = 1; otherwise the kinetic energy part of the lagrangian will 
not remain diagonal under such a transformation. 

Methods for diagonalizing a 2 x 2 complex nonhermitian matrix are known We have 
employed one of them in our analysis in section 5. Here we extend them to arbitrary cases. 
Consider an arbitrary n x n complex matrix Ti. In general Ti has 2n^ parameters. It is 
convenient to separate 7i in to a hermitian and an antihermitian part 

n = M-iT, (66) 

where 

M = M^ = (67) 

r = rt = (68) 

In the study of unstable particle mixing M and F are the mass and decay matrices, respec- 
tively. 

Let us start from simple cases involving real M and F. These correspond to situations in 
which the mixed unstable-particle system conserves CP. Now, Ti is symmetric. It has n(n+l) 
independent parameters. We can diagonalize such a 7i by an orthogonal transformation 



ono 



T 



fx, \ 



V Xn ) 



(69) 



where 00^ = 1. O differs from a usual orthogonal matrix by having complex mixing angles, 
i.e., 9i = Re9i + lm9i {i = l,...,n(n — l)/2). It should be emphasized that a complex 
mixing angle is not equivalent to a real mixing angle times a real phase. Indeed, |cos^| 
could be bigger than 1 if 6' is complex, whereas | cos | < 1 if both 6 and a are real. 
Nevertheless, trigonometric functions with complex angles have very similar properties, such 
as sin^ 9 + cos^ 9 = 1 etc., as the elementary trigonometric functions. This is the advantage 
of introducing complex mixing angles. 

The n{n + 1) parameters of 7i determine the n complex eigenvalues, which have 2n 
parameters, and the n{n — l)/2 complex angles involving n{n — 1) parameters. 

This idea can be extended to situations in which both M and F are complex, provided 
the phases involved in the diagonalization are chosen as complex as well |P0||. 
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Let us first review the diagonalization of a liermitian matrix H 



UHU-^ 



/A, 
V 



\ 

An / 



(70) 



wliere Aj are real. U ^ = U'^ is unitary with n{n — l)/2 real angles and n{n + l)/2 real 
phases. However, not all phases of U are determined by ([70|). To be specific, we write 



U = KU' 



(71) 



where U' is a reduced unitary matrix with n{n — l)/2 angles and n{n — l)/2 phase, K is 
diagonal with n phases 



K 



V 



\ 



(72) 



Clearly, if U diagonalizes H so does its reduced matrix U'. 



An arbitrary Ti has 2n^ parameters: among them, in M and the rest in F. The 



transformation that leaves the unity matrix invariant and diagonalizes Ti. is 



vnv-^ 



V 



\ 

An / 



(73) 



where Aj are complex. In practice we may parameterize V in terms of a reduced unitary 
matrix, such as U' discussed above, but change all the mixing angles and phases into complex 
variables. Now, the 2n^ parameters determine the n complex eigenvalues which involve 2n 
parameters, the n{n — l)/2 complex mixing angles (with n{n — 1) parameters) and the 
n{n — l)/2 complex phases, which also have n{n — 1) parameters. 

For example, for a 2 x 2 7^, has one complex mixing angle and one complex phase (see 
(^)). For a 3 X 3 7Y, has three complex mixing angles and three complex mixing phases. 
It may be parameterized as 
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where c,- = cos 9i and s,- = sin 9j 
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